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ABSTRACT 


The use of the Maximum Principle for the large angle slewing of LSS usually results in 
(he so-called two-point boundary-value problem, in which many requirements (e.g., 
minimum time, small amplitude, and limited control power, etc.) must be satisfied 
simultaneously. The successful solution of this problem depends laigely on the use of 
an efficient numerical algorithm. There are many candidate algorithms available for 
this problem (e.g., quasilinearization, gradient, etc.). Here we discuss only the 
quasllinearization method which has been used for several cases of large angle slewing 
of LSS. The basic idea of this algorithm is to make a series of successive 
approximations of the solution from a particular solvable case (linear or nonlinear) to 
a more general practical case. 

For the rigid spacecraft slewing problem with no constraints on the controls, the 
solution procedure can be found in the literature. This procedure needs to be modified if 
a minimum lime for the slewing problem is desired with control limits given. Recently, 
an indirect method for finding the minimum time was developed to meet all these 
requirements. 

For the general mixed (including both rigid and flexible parts) problem, an additional 
constraint of small vibrational amplitude on the flexible parts is Imposed. To solve 
this problem several steps in which the complexity increases gradually are needed, l.e., 
from a linearized version to a final nonlinear problem, from a less constrained case for 
the control to a more constrained one. from a nonminimum-time level to a near- 
minlmum-tlme slewing in which a trade-off needs to be made between minimum time 
and small flexural amplitude requirements. Some examples of these algorithms are 
presented for planar slewing maneuvers of the SCOLE configuration. 
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Abstract 

The use of the Maximum Principle for the large angle 
slewing of LSS usually results in the so-called two-point 
boundary-value problem, in which many requirements (e.g., 
minimum time, small amplitude, and limited control power, etc) 
must be satisfied simultaneously. The successful solution of 
this problem depends largely on the use of an efficient 
numerical algorithm. There are many candidate algorithms 
available for this problem (e.g., quasilinearization, gradient, 
etc.). Here we discuss only the quasilinearization method which 
has been used for several cases of large angle slewing of LSS. 
The basic idea of this algorithm is to make a series of 
successive approximations of the solution from a particular 
solvable case (linear or nonlinear) to a more general practical 
case. 


For the rigid spacecraft slewing problem with no 
constraints on the controls, the solution procedure can be 
found in the literature. This procedure needs to be modified if 
a minimum time for the slewing problem is desired with control 
limits given. Recently, an indirect method for finding the 
minimum time is developed to meet all these requirements. 

For the general mixed (including both rigid and flexible 
parts) problem, an additional constraint of small vibrational 
amplitude on the flexible parts is imposed. To solve this 
problem several steps in which the complexity increases 
gradually are needed, i.e., from a linearized version to a 
final nonlinear problem, from a less constrained case for the 
control to a more constrained one, from a non-minimum-time 
level to a near-minimum-time slewing in which a trade-off needs 
to be made between minimum time and small flexural amplitude 
requirements. Some examples of these algorithms are presented 
for planar slewing maneuvers of the SCOLE configuration. 
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INTRODUCTION 


MAXIMUM PRINCIPLE IS APPLIED TO 

THE ATTITUDE MANEUVER AND VIBRATION CONTROL 

OF LARGE SPACE STRUCTURES 

(A) PERFORMANCE INDICES 

(B) BOUNDARY CONDITIONS 

(C) CONTROL REQUIREMENTS 

THIS LEADS TO THE TWO-POINT BOUNDARY-VALUE PROBLEM 
(TPBVP) 

ONE OF THE METHODS OF S0LVIN6 TPBVP IS THE 
QUASILINEARIZATION ALGORITHM 
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MAXIMUM PRINCIPLE 


STATE EQUATIONS 

x - f(x) ♦ B(x)u, x(0)-x 0 , x(t f )-x f ( I ) 

PERPORMANCE INDICES 


J,=(l/2) J (x T Qx ♦u T Ru)dt (2) 

0 

J 2 = I OWt* t f lujU u ib , i=l ... n (3) 

0 

NECESSARY CONDITIONS 

H t =(1/2)(x T Qx ♦ u t Ru) ♦ A T (f(x) ♦ Bu) (4) 

A - - (dHj/dx), A (0) unknown (5) 

(3H|/du)=0, Ru=-B T A (6) 

H 2 = I ♦ A T (f(x)*Bu) (7) 

A = - (dH 2 /dx), A (0) unknown (8) 

u i* - Ujjj sign(B T A ) , 1*1 ... n (9) 

TPBVP 

z = g(z) , z=l x , A ] T = [ Z| t z 2 1 T (10) 

Zj(0). Z|(t f ) known; 
z 2 (0), z 2 (tf) unknown. 
z 2 (0) to be determined. 
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QUASILINEARIZATION ALGORITHM 


(A) LINEAR DIFFERENTIAL EQUATION: 

Nonhomogeneous: z = Az ♦ B. z=[zj. z 2 ] T . (11) 

Zj(O), z | (t f ) known, z 2 (0) to be determined 
Homogeneous: z - Az (12) 

(a) n solns. of (12) ♦ 1 particular soln. of (1 1) 

(b) n ♦ 1 particular solns. of ( 1 1 ) 


(B) NONLINEAR CASE: 


Linearized equation of (10): 

z:(k+ * * =(dg/dz) z^* 1 ^ ♦ h( z^ ) (13) 

where 

z^ is the k* h approximate solution 

of the nonlinear equation (10), 
z (k+1) sZ (k) +Az (k) 

Z“[z | , Z 2 ] T , 

Z|^ k+ ^(0). Z| (k+I >(t f ), known 
z 2 (k+1) (o) to be determined 
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Spacecraft Coma. Lab Exp&imsnt 

(SCOLE) 
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TOWARDS 



X Cu 
H O 
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Fig. lb. Attitude of the SCOLE Showing Antenna Line of Sight 











PLANAR SLEWIN6 OF FLEXIBLE SCOI F 


LINEARIZED EQUATION OF MOTION- 


1 -T 


| 1 ^ 

e 

♦ 

0 

a T 

l « 

e 

= 

■ M 

ft 


• t 

i 


o 

c 

7 



I 


0 





<P | <t>2 


U 1 

u 2 



where 

8 Is the angle of rotation, 

7 nxl is the amplitude vector of the flexible modes, 
n is the number of mode used, 

I is the moment of inertia about the axis of rotation 
m. M are the Inertia parameter vector, matrix. 

T Is the stiffness matrix, 

<£(z> is the mode shape function vector. 

0 , Zj is the coordinate along z axis, 

L Is the length of the beam, 
u s is the control torque on the Shuttle, 

Uj are the control actuators on the beam and the 
reflector. 
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STATE EQUATIONS 


s - As ♦ Bu 


s = 


s l 

I S2 


si = 


e 

1 


s2 = 


BOUNDARY CONDITIONS FOR s 


s(0)= 


e f 


0 


0 

. S(t f )= 

0 


0 


0 


0 


0 

2(n* 1 )x 1 


where n is the number of mode shapes used. 
PERFORMANCE INDEX 

J =( I /2) j (x T Ox ♦u T Ru)dt 
0 


'X^S 


TPBVP 


z = Cz . z=l s .A ] T = l Z|. z 2 1 T 

A is the costate vector, 

Z|(0), z | (t j) known; 
z 2 (0) to be determined. 
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NUMERICAL RESULTS 


(A) SLEWING ABOUT X-AXIS (ROLL) 


0„lu q t ( -first mode. sUfe ) IS used. 

&f ~zoV'j), tf - **(•■> 

y (s T QS + U T R «) dt s« ie,%, 7‘ 1 



CASE / : V, ts usacI , 


cASe 2 : 


Uj 


case 3 : 



CASE 3 •' 


at, Uj acre. usee/ 



676 


CETA- (DEG) 
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CASE: U1 CONTROLS USED 



MODE 1 - ( FT ) 
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CASE: U 1 CONTROLS USED 


UI-XIOOOO(FT-LB) 
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CASE: U*( 1 0 0 0 ). Q*( 0.000 0.001 0.000 0.000 



UI-XIOOOO(FT-LB) 






♦ ♦ » » ♦ 

m yj Qj yj yj 

o> n own 

« m oo t** n 
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i i 
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-0. 39E+00 




0 . 746*01 



0 . 746*01 






U1*X10000(FT-LB) 



CASE: U-( 1 2 6 O ). 0-( 0.000 0.010 0.010 0.010 


U2-X 10( LB) 



-0.22E+00 
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UI-XIOOOO(FT-LB) 
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-0.34E+00 








U3-X 10(18) 
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CASE: U*( 10 3 0 ). 0-( 0.000 0.010 0.010 0.010 



CONCLUDING REMARKS 


1) Solution has been obtained for nonlinear rigid spacecraft 
attitude maneuver (including the rigidized SCOLE). 

2) Use of the Maximum Principle can make the states 
satisfy the boundary conditions very well. 

3) Due the fact that the costates must be used in the method, 
the dimension of equations of the system is doubled, and 
higher computational ability is needed in this method. 

4) Further work on more complicated models (nonlinear 
differential equation) is needed. 

5) Need to consider different cost functions and perform 
parametric studies. 
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